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This paper studies a Stieltjes-type moment problem defined by the 
generahzed lognormal distribution, a heavy-tailed distribution with ap- 
plications in economics, finance and related fields. It arises as the dis- 
tribution of the exponential of a random variable following a general- 
ized error distribution, and hence figures prominently in the EGARCH 
model of asset price volatility. Compared to the classical lognormal dis- 
tribution it has an additional shape parameter. It emerges that moment 
(in)determinacy depends on the value of this parameter: for some val- 
ues, the distribution does not have finite moments of all orders, hence the 
moment problem is not of interest in these cases. For other values, the 
distribution has moments of all orders, yet it is moment-indeterminate. 
Finally, a limiting case is supported on a bounded interval, and hence de- 
termined by its moments. For those generalized lognormal distributions 
that are moment-indeterminate Stieltjes classes of moment- equivalent 
distributions are presented. 
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1. INTRODUCTION 

The moment problem asks, for a given distribution with distribution function (CDF) F 
with finite moments mk{F) = x'^ dF{x) of all orders k = 1,2,..., whether or not 

F is uniquely determined by the sequence of these moments. If F is uniquely determined 
by this sequence, F or a random variable X following this distribution are called moment- 
determinate (for brevity, M-det); otherwise F or X are called moment-indeterminate (M- 
indet). Cases where the support of the distribution F is the positive half-axis IR+ = [0, oo) 
are called Stieltjes moment problems, cases where the support is the real line are called 
Hamburger moment problems, and cases where the support is a bounded interval are called 
Hausdorff moment problems. 

The probably most widely known example of an M-indeterminate distribution is the log- 
normal distribution, first described by Stieltjes (1894/1895) in a non-probabihstic setting 
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and further developed by Heyde (1963). The lognormal distribution is a basic model for de- 
scribing size phenomena in economics and related fields (see, e.g., Kleiber and Kotz, 2003), 
including distributions of personal income, actuarial losses, or city sizes. It also arises in 
mathematical finance in the fundamental geometric Brownian motion model of asset price 
dynamics. Given the central role of the lognormal distribution in Stieltjes-type moment 
problems it is, therefore, of special interest to explore closely related distributions with 
respect to M-indeterminacy. Recently, Lin and Stoyanov (2009) studied a generalization of 
the lognormal distribution derived from a skewed generalization of the normal distribution, 
finding that it is M-indeterminate for every value of the skewness parameter. The present 
paper explores a family of generalized lognormal distributions derived from a more clas- 
sical symmetric generalization of the normal distribution, which compared to the normal 
distribution has an additional shape parameter. Like the classical lognormal distribution, 
this generalized version has been employed in financial economics as well as in modeling 
size distributions. 

It turns out that this family of distributions sheds new light on the classical lognormal 
moment problem, in that M-determinacy now depends on the value of the shape parameter. 
Specifically, the family incorporates heavy-tailed distributions for which not all integer 
moments exist, moderately heavy-tailed distributions for which all moments exist yet the 
distributions are M-indeterminate, and, limiting distribution with bounded 

support that is, therefore, determined by its moments. It also emerges that the classical 
lognormal distribution does not constitute an extreme case within the family: in the setting 
considered here, there exist more as well as less heavy-tailed M-indet distributions than 
the lognormal. 

The paper is organized as follows: Section 2 provides some background on the generalized 
lognormal distribution. Section 3 contains a characterization of moment (in)determinacy 
for the family of generalized lognormal distributions in terms of their shape parameter, 
while Section 4 describes Stieltjes classes pertaining to the indeterminate cases. Section 5 
concludes. 

2. THE GENERALIZED LOGNORMAL DISTRIBUTION 

Being one of the basic distributions in probability and statistics, the normal distribution 
has triggered a number of generalizations. One such generalization is defined by the density 

m = 2 r^/r a r(l + 1/r) {"^^"^ " ^"^ ' "^^ < ^ < ^' (2-1) 

which includes the normal as the special case where r = 2. Here yU G M is a location 
parameter and a G M"*" is a scale parameter. The new parameter r G M"^ is a shape 
parameter measuring tail thickness, with lower values of r indicating heavier tails. The 
parameter r plays a crucial role below. 
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This distribution is fairly widely known; however, it is known under different names in 
different fields and it was (re) discovered several times in different contexts. Specifically, 
since r = 2 yields the normal distribution and r = 1 the Laplace distribution, the distribu- 
tion (2.1) is known both as a generalized normal distribution, in particular in the Italian 
language literature (Lunetta, 1963; Vianelli, 1963), and as a generalized Laplace distribu- 
tion. It is also known as the normal distribution of order r, again especially in the Italian 
literature (e.g., Vianelli, 1983), and as the generalized error distribution, notably in econo- 
metrics and finance (e.g.. Nelson, 1991). A further name is exponential power distribution 
(Box and Tiao, 1973), the name under which this distribution is presumably best known 
in the statistical literature. To the best of the author's knowledge, the generalized form 
(2.1) was first proposed in a Russian journal by Subbotin (1923), who sought an axiomatic 
basis for a generalized form of Gauss's "law of error." Hence the name Subbotin distri- 
bution is also in use, notably in econophysics (e.g., Alfarano et al., 2012). A multivariate 
generalization of (2.1) is the Kotz-type distribution (Kotz, 1975). 

In what follows we sometimes set fi = 0, since in the context of moment problems 
no extra generality is gained by including this location parameter. There exist different 
parameterizations of (2.1), notably regarding the scale parameter, but for the purposes of 
this paper the relevant parameter is r, so this complication shall be ignored below. 

The generalized lognormal distribution (Vianelli, 1982, 1983), or perhaps logarithmic gen- 
eralized normal distribution, is less widely known than the generalized normal distribution. 
In fact, most of the currently available works are written in Italian and published in Italian 
journals and collected volumes that are often not easily available outside of Italy. A more 
accessible source may be Kleiber and Kotz (2003, Ch. 4.10), who summarize many basic 
properties. The distribution is defined as the distribution of X = exp(y), where Y follows 
eq. (2.1), leading to the density 



If a random variable X follows eq. (2.2) this is denoted as X ~ GLN(/i, a, r). The 
distribution will sometimes be referred to as the generalized lognormal distribution of 
order r if further emphasis is needed. The case where r = 2 gives the classical lognormal 
distribution. In eq. (2.2), is a scale parameter, while a and r are both shape parameters. 
The effect of the new parameter r is illustrated in Figure 1. This Figure suggests that 
the density becomes more and more concentrated on a bounded interval with increasing 
r. Specifically, for r = 1.5 the density is much like the classical lognormal density, but 
with slightly heavier tails, while for r = 15 several points of inflection and a more rapid 
decrease in the tails emerge. The limiting case where r — )■ oo will also be explored below, 
see Theorem 3. 




< X < oo. 



(2.2) 
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X 



Figure 1. — Some generalized lognormal distributions (solid grey: fi = 0, a = 1, r = 1.5, 
dashed-dotted grey: fi = 0, a = 1, r = 15). The dashed black curve corresponds to the 
classical lognormal distribution (r = 2, with /i = 0,0" = 1). 



Like the classical lognormal distribution, the generalized lognormal distribution has been 
employed in economics and finance. As mentioned above, it has been used as a model for the 
size distribution of personal incomes. In an application to Italian income data, Brunazzo 
and PoUastri (1986) estimate r in the vicinity of 1.45, suggesting a model with even heavier 
tails than the classical lognormal distribution for their data. It will emerge below that their 
estimated model is not determined by its moments. 

Perhaps more prominently, the distribution also arises in the widely used exponential 
GARCH (EGARCH) model of asset return dynamics (Nelson, 1991), where it provides a 
more realistic specification of the innovation distribution in the volatility equation than the 
normal distribution. Recall that, in view of the exponential transformation employed in the 
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EGARCH model, a widely used alternative to the normal distribution in GARCH modeling, 
the t distribution, leads to tails that are too heavy, in the sense that the distribution 
corresponding to the exponentiated random variable has no moments of any order. In 
contrast, it will emerge below that the less extreme members of the generalized lognormal 
distribution possess moments of all orders, yet they are M-indeterminate. Specifically, all 
models estimated by Nelson (1991), with shape parameters r in the vicinity of 1.56-1.57, 
are not determined by their moments. More recent work (e.g., Taylor, 2005) confirms that 

I < r < 2 is the empirically relevant range of the tail thickness parameter in this model. 
All of these objects are M-indeterminate. 

3. GENERALIZED LOGNORMAL DISTRIBUTIONS AND THE MOMENT PROBLEM 

How can one determine whether or not a given distribution with CDF F is determined 
by the sequence of its moments? Although necessary and sufficient conditions are known 
(see, e.g., Shohat and Tamarkin, 1950), they are not very practical. For M-determinacy, a 
sufficient condition is the existence of the moment generating function (MGF) mx{t) = 
E[e*^] = /~ e*^ dFx{x), |t| < to, for some to > 0. 

From the expression for the density (2.2) of the generalized lognormal distribution it is 
immediate that, for any r G M"*", E[e*"^] = oo for all t > 0; hence the MGF does not exist. It 
remains to explore the existence of the moments themselves. (Note that in view of X > 
(a.s.) it is possible to consider moments of fractional order.) Without loss of generality, set 

II = since exp(yu) is a scale parameter. Substituting z = \nx yields, for some C > 0, 

roo roc 

E[X''] = / x^f{x)dx = C exp{A;z- |z|7(m")} dz. (3.1) 

This shows that convergence of the integral depends on the value of r: for r > 1 the 
integral is finite for all k, for r = 1 the condition \k\ < \ jo is needed, while for r < 1 it 
does not converge for any A; 7^ 0. The following proposition collects these observations: 

Proposition 1 Suppose X ~ GLN{n,a,r). 

(a) The moment-generating function of X does not exist for any r G (0, 00). 

(b) The kth moment E,[X^] exists if and only if 

• k = 0, if r < 1. 

• \k\ < 1/a, if r = 1. 

• k E (—00, 00), if r > 1. 

Apart from the integral representation (3.1), it is also possible to obtain series expansions 
of the moments (when they exist). For r > 1, they are of the form 
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{2^)\ 



2i 



r 




2i + 1 



r 



) 



fc = 0,l,2,..., 



see Brunazzo and Pollastri (1986) or Nelson (1991)^. 

In view of Proposition 1 not all generalized lognormal distributions are of interest in 
the context of the moment problem. For r = 1, only some moments exist, for r < 1 no 
moments exist. The cases where r < 1 therefore provide examples of distributions without 
any moments, integer or fractional. An earlier example was given by Kleiber (2000). For 
the remaining cases where 1 < r < cxd all the moments are finite yet the MGF does not 
exist. These are circumstances under which M-indeterminacy may arise. 

It remains to show that the distributions where 1 < r < oo are indeed M-indet. For 
M-indeterminacy, a useful sufficient condition is the Krein condition (e.g. Stoyanov, 2000). 
In a Stieltjes-type moment problem, it requires, for a density / that is strictly positive for 
all x > a > 0, for some a > 0, that the normalized logarithmic integral of the density 



is finite. Ks[f] is called the Krein integral of /. 

The following Theorem shows that generalized lognormal distributions of orders 1 < r < 
oo are M- indeterminate: 

Theorem 2 All generalized lognormal distributions GLN{fi, a, r) of order 1 < r < oo 
are M-indeterminate. 

Proof. Setting without loss of generality ^ = and a = 1, the Krein integral (3.2) is, 
for a > and Cr > the normalizing constant. 



Since for large x the integrand is eventually dominated by x ^ for any 6 G (0, 1), this 



Alternative proofs could employ results presented by Gut (2002, Remark 6.2) or Pakes 
et al. (2001, p. 110). 



^It should be noted that these works employ different parameterizations of the distribution. Also, Nelson 
(1991) obtains expectations of somewhat more general objects. Setting 7 = 0, p — and 6* = 1 in his 
Theorem A1.2 yields the required moments. The resulting expressions can be shown to coincide with those 
presented by Brunazzo and Pollastri (1986). 




(3.2) 




integral is finite for all 1 < r < 00, which gives the result. 



□ 
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For Xi ~ GLN(/ij, (Tj, Tj), i = 1,2, with > 1 and densities /j it is easily seen that 
hm;j._j.oo fi{x)/f2{x) = oo iff ri < r2, hence the generahzed lognormal distributions are, in 
a sense, "more M-indeterminate" for smaller r. (Indeed, in view of Proposition 1 for r = 1 
some moments no longer exist.) Specifically, the generalized lognormal distributions with 
1 < r < 2 are even more extreme than the classical lognormal distribution (r = 2). Also, 
the cases where 2 < r < oo are less extreme. It is also worth noting that although the tails 
of the generalized lognormal distribution become lighter and lighter with increasing r, the 
distribution is M-indet no matter how large r. It is, therefore, natural to ask what happens 
in the limit, i.e., for r — )■ oo. The following Theorem addresses this case: 

Theorem 3 For r — )■ oo, the generalized lognormal distribution GLN(fi,a,r) tends to a 
distribution supported on a bounded interval. Hence this limiting distribution is M-det. 

Proof. It is convenient to analyze the limiting case for the distribution of Y = InX, 
i.e., the generalized normal distribution. Without loss of generality, set fi = and a = 
1. A random variable Y following a generalized normal distribution admits the mixture 
representation (Devroye, 1986, p. 175) 



Y = UZ (3.3) 

where U is uniform on [—1, 1] and Z ~ (r y^)W'^/'- with W ~ Ga(l + 1/r, 1), i.e, a gamma 
distribution with scale 1 and shape parameter 1 + 1/r. Hence Z follows a generalized gamma 
(GG) distribution, specifically Z ~ GG(r, r^/^', 1 + 1/r). The moments of Z are (see, e.g., 
Kleiber and Kotz, 2003, p. 151) 



^.^k. _ (rV)WAr(l + (fc + l)/r) 

^z \ - r(i + i/r) ' ^-i'^'---- 

Now limj._s.oo IE[Z'^] = 1 for all k, and it follows that Z = r^^'^W^^^ tends to a point mass 
at 1 by Frechet-Shohat (e.g., Galambos, 1995, p. 81). Thus limr_s.oo Y = U, and the density 
of exp(U) is given by 



fix) = e-'<x< e. (3.4) 

2 X 

This distribution has compact support, hence it is determined by its moments. □ 



Lunetta (1963) provides an alternative derivation of the limiting distribution of the gener- 
alized normal distribution that analyzes the limit of its characteristic function. However, we 
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prefer the approach involving a mixture representation presented here because it motivates 
further questions, on which more below. 

Interestingly, Bomsdorf (1977) observed that a distribution of the type described by 
eq. (3.4) occurs as the distribution of prizes in lotteries, hence he calls it the prize compe- 
tition distribution. Among other characteristics he also provides the MGF of this object. 

4. STIELTJES CLASSES FOR MOMENT-INDETERMINATE 
GENERALIZED LOGNORMAL DISTRIBUTIONS 

The preceding section showed that generalized lognormal distributions of orders 1 < 
r < oo are M-indeterminate, by way of an existence proof. To round off the discussion, 
this section provides explicit examples of distributions that are equivalent, in the sense of 
having identical moments of all orders, to these indeterminate distributions. 

A Stieltjes class - a term coined by Stoyanov (2004) - corresponding to a moment- 
indeterminate distribution F with density / is a set 



S{f,p) = {fs{x)\fe{x):=f{x)[l + ep{x)],xesnpp{f),6e[~l,l]}, 

where p{x) is a perturbation function satisfying —1 < p{x) < 1 and E[X'^j9(X)] = for all 
fc = 0,1,2,.... 

It is possible to obtain Stieltjes classes for the generalized lognormal distributions of 
orders 1 < r < oo that generalize a recently derived Stieltjes class pertaining to the 
classical lognormal distribution. The construction of the required Stieltjes classes in the 
following Theorem is adapted from a construction presented by Stoyanov and Tolmatz 
(2005, Theorem 3): 

Theorem 4 Suppose X ~ GLN{ix, a, r) with density fr, (/i, a, r) G M x IR+ x (1, oo). 

(a) The function 

J sin{(x - l)^/^}exp {^1 Inx - + Inx - (x - 1)^/^} , a; > 1, 
^^""^ ~ 1 0, a;<l,^ ^ 

is bounded on IR+ for all {n.a.r) e M x M+ x (l,oo), with ¥.[X^hr{X)] = for all 
fc = 0,1,2,.... 

(b) pr '■= hr/Hr, with Hj. := sup^ |/ir-(a^)|; defines a perturbation corresponding to fr- 

(c) The family of functions fr,e{x) = fr{x)[l + e Pr{x)], e G [—1,1], defines a Stieltjes 
class comprising distributions whose moments are identical to those of for any 
ee [-1,1]. 

Proof. The function hr is continuous on (1, oo), with lima;_j.i+ hr{x) < oo and lim^-^oo hr{x) 
0, hence is bounded on M"*". 
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By construction, with > the normahzing constant of fr, 



for = 0, 1, 2, . . . in view of Lemma 1 of Stoyanov and Tolmatz (2005) and the fact that 



This proves (a). 

Since Hr := sup^ |/?-r(a;)| < oo we may set Pr{x) = hr{x)/Hr, assuring |pr(a;)| < 1 for all 
X. This gives (b). Finally, fr,e{x) = fr{x)[l + e Pr{x)] defines a density for any e G [—1,1], 



It should be noted that the construction of Stoyanov and Tolmatz (2005) is somewhat 
more general, in that the kernel k{x) := {x — 1)^/^ used here may be generalized to 
a three-parameter family of kernels defined by k{x]^,6,/3) := {6x — ^)^ tan(7r/3), where 
(^, 5, /3) G X ]R+ X (0, 1/2). Thus amending the kernel in this manner defines a four- 
parameter family of perturbations pr{x;^,S, (3) leading to Stieltjes classes that generalize 
the three-parameter family of Stieltjes classes for the classical lognormal distribution de- 
rived by Stoyanov and Tolmatz (2005). However, the Stieltjes class presented above already 
provides infinitely many distributions whose moments coincide with those of the generalized 
lognormal distribution. 

In (4.1), the choice of f3 = 1/4 was made because it is related to one of the classical 
examples of an M-indeterminate distribution that dates back to the pioneering work of 
Stieltjes (1894/1895). Stieltjes considered the case where ^ = and the perturbation h{x) = 
sin(a;^/'^), x > 0, used in the proof of part (a) of Theorem 4; it pertains to a certain 
generalized gamma distribution. Moreover, a shift ,^ > is needed in (4.1), as otherwise 
the resulting object would exhibit a singularity at the origin, see also the discussion in 
Stoyanov and Tolmatz (2005, Section 4). 



The paper exhibited a family of distributions, occurring in economics and finance, that 
generalizes the lognormal distribution, the classical example of a moment-indeterminate 
distribution. It emerged that a large subfamily consists of moment-indeterminate distri- 
butions, but also that not all members share this property of the lognormal, for different 





which gives (c). 



□ 



5. FURTHER DISCUSSION AND CONCLUDING REMARKS 
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reasons: some tails are so heavy that not enough moments exist, while a limiting case 
corresponds to a light-tailed distribution with compact support. 

It may, therefore, be asked to what extent it is possible to characterize the gener- 
alized lognormal distributions with r = 1, i.e. the log-Laplace distributions, for which 
E[X'^] < oo iff |A;| < l/a. If one leaves the classical setting of the moment problem 
characterizations in terms of certain moments are possible. First, Th. 1 of Lin (1992) im- 
plies that characterizations in terms of fractional moments are feasible: for a sequence 
{kn I < kn < l/a;n G N} of positive and distinct numbers converging to some 
kf) G (0, l/o"), the sequence {E[X'^"] | n G N} of fractional moments characterizes the 
distribution. Second, observe that for r = 1 the first moment exists iff cr < 1. It is well 
known that existence of the first moment permits characterization of the underlying distri- 
bution in terms of the triangular array of first moments of the associated order statistics, 
{E[Xfc.„] I k = 1,2, ... ,n;n G N}, where Xi,n < X2;n < • • • < -^n:n are the order statis- 
tics in a sample of size n. In fact, certain subsets of this array are already sufficient, see 
Huang (1989) for a review. Such characterizations are meaningful in applications to income 
distribution (Kleiber and Kotz, 2002), one of the fields where the generalized lognormal 
distribution has been employed. Note also that both characterizations, via fractional mo- 
ments as well as via moments of order statistics, are available for all generalized lognormal 
distributions with r > 1 since moments of arbitrary order exist in that case. 

It is natural to ask about M-determinacy of the more widely known distribution of InX, 
the generalized error or Subbotin distribution (2.1). This is a Hamburger moment problem. 
The answer is already available in the literature, although not in a probabilistic setting: the 
family of generalized error distributions also admits M-indet examples, namely for r < 1, 
and a Stieltjes class is given in Shohat and Tamarkin (1950, p. 22). 

It is also known that for some M-determinate distributions power transformations lead to 
M-indeterminacy and vice versa (e.g. Stoyanov, 1997). The standard example is the general- 
ized gamma distribution. For X ~ GLN(/i, a, r), it is easily seen that ~ GLN(p/i,p(T, r) 
for all p > 0, showing that the distribution is closed under power transformations. Hence 
this well-known property of the classical lognormal distribution extends to the generalized 
version (2.2). Consequently, consideration of power transformations does not lead to new 
insights regarding the moment problem here. 

However, it might be worthwhile to further explore aspects of the mixture representation 
(3.3). This representation is a special case of a general mixture representation for uni- 
modal distributions known as Khinchine's theorem. The exponentiated version states that 
exp(y) = exp(?7Z), i.e. a random variable following a generalized lognormal distribution 
can be obtained as the exponential of the product of a uniform and a transformed gamma 
random variable. It would be interesting to characterize the set of mixing distributions Fz 
leading to indeterminate log-unimodal distributions. 
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